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Abstract 



Multiplicities of charged secondary hadrons in the relativistic nucleus- 
nucleus collisions in a wide energy range are analysed on the basis of partial 
stimulating emission and cluster cascading models. The experimental data 
are obtained by means of the two metre propane bubble chamber of JINR 
(Dubna). The results are compared with the corresponding data for pp and 
pp-collisions at higher energies. It is shown that the regularities in pp and 
pp-interactions at high energies and wide energy range and in nucleus-nucleus 
interactions at relatively low energies per nucleon and narrower energy range 
are similar. 

Key words: relativistic heavy ions, multiplicity distribution, partial stimu- 
lating, cluster cascading. 
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1 Introduction and Basic Notions 



The interest to the study of multiplicity distributions of charged secondaries 
has been increased again after the new powerful accelerators have been con- 
structed and the beams of protons, antiprotons, electrons, positrons and 
heavy ions have been obtained (see,e.g. [1-5]). 

In the present work multiplicities of charged secondary hadrons in pp, pp 
and nucleus-nucleus collisions in a wide energy range are analysed on the basis 
of partial stimulating emission (PSE) and cluster cascading models (CCM). 
The manifestation of the negative binomial distribution (Pascal distribution) 

Pn = (j^V(^L T Y(^) t (1) 



(A;-l)!n! \<n>+kj \< n > +k 

in the multiplicity distributions of secondaries is interpreted in the framework 
of these models. 

In Eq.(l) < n >=< n± > is the average multiplicity of all charged 
secondaries. Parameter k determines the form of the distribution, e.g. if 
k~ x = 1 we get the geometrical distribution, if k^ 1 = we get the Poisson 
distribution. 

It is useful to establish a recurrance relation between P n and P n+1 [6]. 
When deriving this relation one starts from the consideration that the event 
with multiplicity (n + 1) can be expressed by means of (n + 1) number of 
events with multiplicity n: 

s(n) = ( !i± i^ ±1 (2) 

Inserting Eq.(l) into Eq.(2) one can express g{n) as: 

g{n) = a + bn, (3) 

where: 

k < n > < n > 
a = , b = (4) 

k+ < n > k+ <n> 



From here one can write: 



a , a 

< „ > = — * ( 5 ) 
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Dispersion D = (< n 2 > — < n >2) 1 ^ 2 of the distribution (1) and 
parameters < n > and k are related as: 



k = 



<n> 2 



(6) 



D 2 - <n> 



PSE model admits the following interpretation of the distribution (1) and 
the relation (3) [6]. The emitted particles are uniformly distributed among 
k cells. These cells do not correlate and there is no connection between the 
particles in different cells. The additional (n + l)-th particle can be emitted 
in the initial act of the collision independently of the already existing n 
particles. This is reflected by the constant term a in the function g(n). 
In this case b = (i.e. k" 1 = 0) and the classical Poisson distribution is 
obtained. But such an emission can be intensified as a result of quantum 
interference effects. The average effect of this intensification is expressed by 
adding the linear term bn {b = const) in the function g{n). It follows from 
the Eqs. (3) and (5) that: 



One can conclude that n/k is the average number of particles among 
already existing n particles which promote the creation of a new (n + l)-th 
particle. So k^ 1 is the relative average fraction of particles which stimulate 
such creation. 

On the other hand the multiple production can be interpreted in the 
framework of CCM. One assumes here that after the collision of high energy 
particles (leptons, hadrons, nuclei) some exited n-particle system is produced 
which is formed as an iV-cluster state. Each of these clusters is formed by 
the particles which are produced directly or indirectly from one particle pro- 
duced in the initial act of the collision. These latter particles are called the 
"patriarchs" of the clusters. The "patriarch" which does not produce secon- 
daries forms one particle cluster itself. It is assumed that "patriarchs" and 
consequently clusters are produced independently from each other. Therefore 
for the multiplicity of clusters the Poisson distribution holds: 




(7) 



F(N) ~ — <N> 



N 



(8) 
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For the better understanding of the CCM let us give a short review of the 
formulae of this model. It is evident that the average number of clusters 
< N > is given by the formula: 

<iV>=^- (9) 

<n c > 

where < n > is the average number of charged hadrons, < n c > is the average 
number of hadrons in the cluster. 

Let F c (n c ) be the distribution of particles inside one cluster. It is assumed 
in the CCM that the recurrence relation between F c {n c + 1) and F c (n c ) is of 
the form: 

(n c + l)F c (n c + 1) 
g c [n c ) = — — - = bn c , b = const, n > 1 (10) 

The meaning of this relation is that the effect of the creation of (n c + 1)- 
th particle is proportional to the number of particles already existing in the 
cluster in average. Iterating Eq.(10) one gets: 

F c (n c )=F c (l) (11) 

n c 

where F c (l) can be found from the normalization condition: 

S>c(0 = l (12) 



W = (13) 



From here one gets: 

One can find the following average values: 

<n c >=J2n c F c (n c ) = ^l (14) 

<n c 2 >=E^eK) = ^f (15) 

Tic 

Note finally that on the basis of Eqs.(8) and (11) one can write the ex- 
pression for the total multiplicity n in the form of the negative binomial 
distribution which is represented in the form [6]: 

P n ~ a(a + b)[a + b(n - 1)] (16) 
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where 



a =< N > F C {1) 



(17) 



Inserting Eq.(16) into the recurrence relation (2) one gets again Eq.(3). 
The identification of the parameters a and b which is given by Eqs.(lO) and 
(17) completly corresponds to their earlier physical meaning. So the CCM is 
compatible with the negative binomial distribution. 

2 Analysis of the Experimental Data 

The experimental data are obtained on the two-metre propane bubble cham- 
ber (PBC-500) of the Laboratory of High Energies of JINR (Dubna) with 
tantal targets in it which were bombarded by p, d, He and C beams[5]. The 
data on pp (ISR and others) and pp (UA5- Collaboration) collisions in a 
wide energy range [1-4,6] are used for comparision. 

It was shown in Ref.[4] that in pp and pp- collisions (10 < y/s < 900Ge^) 
k~ x increases linearly and reaches the value 0.31 at y/s = 900GeV, i.e. every 
third particle is activly interacting. The average number of charged parti- 
cles in the cluster < n c > reaches the value 4.55, the average number of 
clusters reaches the value 8. May be there is even a decrease of the num- 
ber of clusters as compared to the energy y/s = 540GeV (Table 1). The 
width of the distribution of charged particles in the clusters increases with 
increasing energy much faster than the corresponding dependence for the to- 
tal multiplicity. The ratio R =< n > / D 2 decreases 12 times in the range 
8.33 < y/s < OOOGeV. The same quantity for the clusters R c =< n c > / D 2 
decreases 400 times (Tables la and lb, Fig. 1). Such a decrease is observed 
up to y/s = 200GeV and after that R and R c behave almost in the same 
way. After y/s = 200GeV a tendency of the narrowing of the multiplicity 
distributions in the cluster and the total multiplicity is observed. 

What is the behaviour of R = f(E) and R c = f(E) for £>Ta-collisions in 
the range 2-10 GeV? A fast decrease of R and especially of R c is observed up 
to the energy 5 GeV, i.e. the dispersion increases fast (Table 2, Fig. 2) At the 
energies higher than 5 GeV a change of the regime is observed and probably 
they behave as constants. The constant regime in pp and pp-collisions is 
observed from y/s = 200GeV. 

Consider the same dependences for dTa, HeTa and CTa-collisions (Ta- 
bles 3,4,5). For dTa-collisions R C (E) decreases faster than R(E). The situ- 
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ation is the same as in £>Ta-collisions in the same energy range. In the case 
of HeTa and CTa-collisions the behaviours of R C (E) and R(E) are similar. 

It is interesting to note that a similar picture arises, if we consider the 
dependence of R and R c on the atomic weight Aj, of the projectile nucleus. 
For pTa, dTa, HeTa and CTa-collisions at 4.2 AGeV R and R c decrease with 
increasing Aj. But, if R decreases 19 times, R c decreases 50 times. Further, 
the ratio R c j R is approximately equal to 8 for pTa-collisions and is equal to 
1.5 for CTa-collisions at the same energy, i.e. this ratio tends to one with 
incresasing Ai. 

In Ref.[7] the dependence k^ 1 = f(s) for pp and ^-collisions in the range 
10 < \/s < 900GeV is approximated by the formula: 



A;" 1 > and decreases rather slowly with increasing energy (bi ~ 0.06). 
It has been shown by our analysis that at lower energies, in particular, at 
y/s < 8.33GeV (corresponding laboratory energy is 36 GeV) k' 1 < 0. At 
higher energies k~ l becomes positive. 

Consider now nucleus-nucleus collisions. If by analogy with Eq.(18) we 
approximate A; -1 by the formula: 



it turnes out that in a narrow energy range (2-5 GeV) a fast increase of k^ 1 
is observed, b 2 > 0.15 (Tables 3,4,5, Fig. 3). For the energies higher than 5 
GeV the fast increase of k~ l is slowed down. The essential feature of pTa- 
collisions is that at the energy 2.48 GeV the parameter k~ x becomes negative. 
Starting from 4.3 GeV (may be earlier) k~ x becomes positive, i.e. there is 
a similarity with pp- collisions but at lower energies. For dTa, HeTa, CTa- 
collisions parameter A;" 1 is positive and increaseing with increasing energy, 
i.e. the growth of the atomic weight of the incoming nucleus plays the same 
role as the growth of energy in pp-collisions (Tables 1, 3-5,Refs.[l-8]). 

Consider the situation with negative A; -1 in some detail. Note that one 
writes nothing about them in Ref. [6], though they are obtained from the 
experimental data (Tables la, 2, Refs.[9,10]). Probably the reason is that the 
negative bimomial distribution and PSE and CCM in their strict formulations 
assume the positivness of k^ 1 . But the stable presence of negative k~ x made 
us to take a broader view to this problem. First of all, if we consider the 



k 1 = ai + b\ In \fs , 



(18) 



k 1 = a 2 + b 2 In E, 



(19) 
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negative binomial distribution in the form of Eq.(16) and compare it with 
the Polya-Egenberger distribution (see, e.g. [10]) 

where g 2 = (D 2 — < n >)/ < n > 2 ,it is evident that these two distributions 
coincide. The role of the parameter k~ l = b/a is played by g 2 . But the sign 
of the parameter g 2 in the Eq.(20) is not fixed. In particular, if g 2 = 0, we 
get the Poisson distribution, if g 2 > 0, this distribution is broader than the 
Poisson one, if g 2 < 0, it is narrower than the Poisson one. So parameter 
g 2 can be thought as a measure of deviation of the distribution from the 
Poisson one and if we perform the analysis of the experimental data in terms 
of the Polya-Egenberger distribution, parameter k^ 1 = g 2 can take positive 
and negative values as well. Further, there are some grounds to think that 
the negative values of k^ 1 are compatible with the assumptions of PSE and 
CCM. In fact, the second term in Eq.(3) corresponds to quantum-mechanical 
interference effects. If this effects stimulate the creation of (n+l)-th particle, 
then b > and hence k~ x = b/a > 0, since a > 0. But it is natural to suppose 
that this effects can make weaker the emission (cupture of the particle). In 
this case j3 < and hence k~ x < 0. So one can conclude that k~ x < 
corresponds to the cupture of secondary particles (see also [5,11]). 

Proceed now to consider the energy dependence of the total average mul- 
tiplicity of particles < n > and the average multiplicity of particles in the 
clusters < n c > in AjTa-collisions and compare this data with correspond- 
ing results on pp and pp-collisions. In the energy range 7.42-200 GeV in 
pp and pp- collisions < n > and < n c > increase rather fast, but < n > 
increases faster than < n c >. In the range (200-900) GeV the increase of 
these quantuties is slower and both of them increase 1.6 times (Table 1). 

Consider now pTa, dTa, HeTa and CTa-collisions. In pTa- collisions the 
energy range can also be divided into two parts. The first interval (2-5) GeV 
and the second one (5-9) GeV. In the first interval we observe a more rapid 
increase of the average multiplicities < n > and < n c > than in the second 
one (Table 3, Fig. 4). So we have qualitativly the same situation as in pp and 
pp- collisions, but the energy range is more narrow and low. In dTa, HeTa 
and CTa-collisions in the energy range (2-5) GeV the average multiplicities 
< n > and < n c > increase similarly with increasing energy. 
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It is interesting to trace the dependence of < n > and < n c > on the 
atomic weight Ai of the incoming nucleus at fixed energy. It is seen from 
Tables 3,4,5 and Fig. 4 that at 2.5 GeV < n c > increases 3 times and at 4.3 
GeV - 4 times. The behaviour of < n > is similar. So the increase of the 
atomic weight of the incoming nucleus in AjA r collisions plays a similar role 
as the increase of the energy. As a result in AjA r collisions the same effect is 
achieved at relativly low energies as compared to hadron- hadron collisions. 
This can be confirmed by the following example: in CTa-collisions at 2.48 
AGeV < n c >= 2.66 and at 4.30 AGeV < n c >= 4.55. The same number 
of particles is contained in the clusters in pp-collisions at = 200GeV and 
y/s = 900GeV, respectively. 

It is interesting to study the dependences of the dispersions on < n > 
and < n c >: 

D = f(<n>), (21) 
D c = /(< n c >) (22) 

It is seen from Fig. 5 that at the energies higher than 5 GeV the data for 
the dependence (21) for hadron-hadron and nucleus- nucleus collisions lie on 
two different curves (curves 1 and 2), but the data for the dependence (22) 
lie on the same curve (curve 3). 

Let us show in conclusion the multiplicity distribution of charged secon- 
daries in dTa-collisions at 5.18 AGeV and its fit according to Eq.(l) (Fig. 
6) with only normalization parameter,which is approximately equal to one, 
X 2 /N = 1.15, N is the number of experimental points. 

The authors are indebted to the staff of two metre propane bubble cham- 
ber of JINR (Dubna) for supplying the data. They would like to thank Ya. 
Darbaidze, E. Khmaladze, N. Kostanashvili, P. Pras, G. Roche, T. Topuria, 
M. Topuridze for interesting discussions. On of the authors (V.R.G.) ex- 
presses his deep gratitude to Bernard Michel and Guy Roche for the warm 
hospitality at the Laboratoire de Physique Corpusculaire, Universite Blaise 
Pascal, Clermont-Ferrand, to V. Kadyshesky, T. Kopaleishvili, H. Leutwyler, 
W. Riihl for supporting his stay at the L.P.C. and to NATO for supporting 
this work. 
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Table 1a 

Characteristics of the Multiplicity Distributions of Charged Hadrons 

in ^collisions 





Energy in the c.m.s. y/s, GeV 


7.42 


8.33 


9.78 


30.4 


62.2 


< n > 


4.56±0.04 


4.78±0.03 


5.32±0.13 


10.7 


13.6 


D 


2.09±0.04 


2.23±0.03 


2.58±0.05 


4.59 


6.01 


<n> /D 2 


1.04±0.02 


0.96±0.04 


0.80±0.03 


0.51 


0.38 


k- 1 


-0.009 


0.008 


0.05 


0.09 


0.12 


b 


-0.044 


0.039 


0.20 


0.49 


0.62 


<n c > 


0.98 


1.02 


1.12 


1.43 


1.69 


D c 




0.14 


0.38 


0.87 


1.26 


<n c > /D 2 C 




51 


7.76 


1.88 


1.06 


<n> /k 


-0.04 


0.04 


0.25 


0.97 


1.66 


<N> 


4.66 


4.69 


4.75 


7.48 


8.05 



Table 1b 

Characteristics of the Multiplicity Distributions of Charged Hadrons 

in pp-collisions 





Energy in the c.m.s. y/s, GeV 


200 


546 


900 


< n > 


21.3±0.80 


29.1±0.9 


34.6±1.2 


D 


10.9±0.4 


16.3±0.4 


20.7±0.6 


<n> /D 2 


0.18±0.02 


0.11±0.01 


0.08±0.01 


k' 1 


0.22 


0.27 


0.31 


b 


0.82 


0.88 


0.92 


<n c > 


2.65 


3.46 


4.55 


D c 


2.77 


4.11 


6.01 


<n c > /D 2 C 


0.34 


0.20 


0.12 


<n> jk 


4.65 


7.88 


10.31 


<N> 


8.07 


8.41 


7.60 
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Table 2 

Characteristics of the Multiplicity Distributions of Charged Hadrons 

in pTa-collisions 





Energy per nucleon in the lab. system, GeV 


2.48 


4.30 


5.48 


9.94 


< n > 


2.98±0.06 


4.73±0.07 


6.12± 0.09 


7.84±0.13 


<n> /D' 2 


1.38±0.05 


0.74±0.05 


0.53±0.05 


0.30±0.03 


k- 1 


-0.08±0.01 


0.07±0.01 


0.15±0.01 


0.31±0.02 


b 


-0.30 


0.26 


0.47 


0.71 


<n c > 


0.86±0.04 


1.18±0.05 


1.41±0.09 


1.97±0.12 


<n> jk 


-0.27 


0.34 


0.91 


2.40 


D c 




0.45 


0.82 


1.71 


<n c > /D 2 C 




5.83 


2.08 


1.47 


<N> 


3.47±0.18 


4.01±0.21 


4.34±0.22 


4.22±0.23 



Table 3 

Characteristics of the Multiplicity Distributions of Charged Hadrons 

in <iTa-collisions 





Energy per nucleon in the lab. system, GeV 


2.48 


4.30 


5.18 


< n > 


4.46±0.07 


7.15±0.10 


7.68±0.21 


<n> /D 2 


0.68±0.03 


0.40±0.05 


0.31±0.03 




0.09±0.01 


0.22±0.01 


0.29±0.03 


b 


0.29 


0.61 


0.69 


<n c > 


1.20±0.06 


1.65±0.09 


1.90±0.12 


<n> jk 


0.40 


1.57 


2.23 


D c 


0.50 


1.23 


1.59 


<n c > /Dl 


4.80 


1.09 


0.75 


<N> 


3.72±0.20 


4.22±0.22 


4.04±0.25 
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Table 4 

Characteristics of the Multiplicity Distributions of Charged Hadrons 

in if eTa-collisions 





Energy per nucleon in the lab. system, GeV 


2.48 


4.30 


5.18 


< n > 


7.64±0.14 


10.81±0.15 


11.79±0.25 


< n > /D' 2 


0.42±0.03 


0.21±0.02 


0.16±0.02 


k- 1 


0.22±0.02 


0.35±0.03 


0.45±0.05 


b 


0.61 


0.79 


0.84 


<n c > 


1.67±0.08 


2.41±0.10 


2.88±0.16 


<n> /k 


1.56 


3.76 


5.33 


D c 


1.22 


2.38 


3.12 


<n c > /D 2 C 


0.89 


0.42 


0.29 


<N> 


4.19±0.14 


4.48±0.23 


4.09±0.22 



Table 5 

Characteristics of the Multiplicity Distributions of Charged Hadrons 

in CTa-collisions 





Energy per nucleon in the lab. system, GeV 


2.48 


4.30 


< n > 


12.73±0.61 


19.75±0.39 


<n> /D 2 


0.18±0.02 


0.08±0.01 


k- 1 


0.36±0.05 


0.56±0.04 


b 


0.82 


0.92 


<n c > 


2.66±0.10 


4.55±0.12 


<n> /k 


4.66 


11.03 


D c 


2.78 


6 


<n c > /D 2 C 


0.35 


0.12 


<N> 


4.77±0.22 


4.29±0.24 
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4 Figure Captions 

Fig.l Energy dependence of R and R c in pp-and pp-collisions. + for R, • 
for R c . 

Fig. 2 Energy dependence of R in pTa(»),dTa(^ r ),HeTa(W),CTa(£)- 

collisions and R c in pTa(o), dTa(y), HeTa(D), CTa(A)-coHisions. 

Fig. 3 Energy dependence of the parameter k~ x in A(Ta- collisions. pTa(«), 
dTa(o), HeTa(+), CTa (y). 

Fig. 4 Energy dependence of < n > in A,{Ta- collisions, pTa(+), dTa (A), 
HeTa(U), CTa(y). 

Energy dependence of < n c > in A{Ta- collisions, pTa(«), dTa(A), 
HeTa(u), CTa(y). 

Fig. 5 Dependence D(< n >) in e + e~, pp, pp, A(Ta- collisions, e + e _ (x), 
pp, pp (•), pTa(A), dTa (y), HeTa (y), CTa (A), curves 1 and 2. 

Dependence D c (< n c >), in e + e~, pp, pp, AjTa-collisions, e + e~(o), pp, 
pp(+), pTa(«), dTa(D), HeTa(W), CTa(d]), curve 3. 

Fig. 6 Total multiplicity distribution of charged secondaries in dTa- colli- 
sions at 5.18 AGeV, the curve corresponds to Eq. (1). 
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Figure 6: 
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